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New experimental data on the partitioning of siderophile and chalcophile elements among metallic 
and silicate phases may be used to constrain hypotheses of core formation in the Earth. Three current 
hypotheses can explain gross features of mantle geochemistry , but none predicts siderophile and 
chalcophile element abundances to within a factor of two of observed values. Either our understanding 
of metal- silicate interactions and j or our understanding of the early Earth requires revision. 


Core formation has been described as perhaps the single most 
important thermal event in the history of the Earth 1 . This 
description is certainly true if the core formed by segregation 
of metal in an initially homogeneous body the size of the present 
Earth, as the heat liberated in such an event would raise the 
temperature of the Earth by —2,000 °C (ref. 2). However, less 
extreme models in which core formation is concomitant with 
accretion appear to be more probable 3 . Regardless of whether 
core formation resulted in catastrophic changes in the Earth, 
better knowledge of the core-forming event would certainly 
elucidate important aspects of the histories of the Earth and 
Moon 4 . 

Geochemists have traditionally divided the chemical elements 
into four classes: (1) lithophile — those elements which prefer 
to be combined with oxygen and which typically reside in rocky 
material; (2) siderophile — elements which prefer to be in the 
metallic state; (3) chalcophile — elements preferring to be associ- 
ated with sulphur; and (4) atmophile — elements so volatile that 
they are typically concentrated in a planet’s atmosphere. 
Obviously, membership in one of these classes is not invariant, 
but is modulated by the ambient conditions. For example, as 
the temperature increases, elements which normally exist as 
solids or liquids will vaporize and become gases, and, in fact, 
the classes of lithophile, siderophile and chalcophile are typi- 
cally subdivided into groups of elements of similar volatility. 

In principle, the record of core formation is contained in the 
siderophile and chalcophile element contents of rocks derived 
from the Earth’s mantle. Certainly these elements are much 
more depleted in the upper mantle than lithophile elements of 
similar volatility, implying that they record a depletion event 
unseen by the lithophile elements — presumably core formation. 
In practice, these elemental abundances have been attributed 
to either gross disequilibrium between the core and upper 
mantle, or mixing of siderophile- and chalcophile-element-rich 
material into the upper mantle after core formation, or both 5,6 . 
In this view, little information about core formation may be 
obtained from mantle siderophile and chalcophile elements 
because their abundances were established either after core 
formation ceased or as the results of processes which never 
approached equilibrium. The evidence cited for disequilibrium 
between mantle and core during core formation is simply that 
the concentrations of siderophile and chalcophile elements in 
mantle-derived rocks are orders of magnitude higher than would 
be expected if these rocks had ever equilibrated with metal 6 . 
Furthermore, many of the siderophile elements in the Earth’s 
mantle exist in chondritic proportions (that is, with elemental 
abundance ratios characteristic of primitive meteorites (chon- 
drites) which have not undergone chemical processing). Thus, 
successful models of core formation and for siderophile and 
chalcophile clement abundances in the Earth’s upper mantle 
must be able to explain both high abundances and the apparently 
unfractionated ratios of some (but not all) element pairs. 


Models of core formation are further complicated by our lack 
of knowledge of the chemical composition of the core. All 
models of the Earth’s core require the presence of a light 
element’ to reduce the overall density 3 . If the light element is 
sulphur, silicon or carbon, then these elements were probably 
dissolved into the metal (which eventually formed the core) at 
low pressure as the Earth grew by accretion, and may be relevant 
to understanding the core-forming process. If the light element 
is oxygen or a mixture of oxygen and lithophile elements which 
have followed oxygen into the core, then the light element must 
have been incorporated at high pressure 4 and may not have 
played a role in establishing the siderophile and chalcophile 
element abundances of the upper mantle. 

The behaviours of siderophile trace elements in complex, 
natural metal-silicate systems are poorly known. This lack of 
knowledge has led us to determine partition coefficients (Ds; 
for definition, see below) between solid metal, sulphur-bearing 
metallic liquid and silicate liquid for a suite of siderophile 
and/or chalcophile elements under controlled laboratory condi- 
tions (T= 1,250-1,270 °C; P^lbar). These partition 
coefficients, together with estimates of the mantle abundances 
of the same elements, may be used to evaluate three physically 
plausible but very different hypotheses of core formation. These 
hypotheses are: (1) inefficient core formation, involving equili- 
brium between solid and liquid metals and silicates, with a small 
fraction of solid and liquid metal remaining trapped in the 
mantle, subsequently to be oxidized 7 ; (2) equilibrium between 
an Fe-S-0 metallic liquid and the mantle 8 ; and (3) 
heterogeneous accretion involving a late ‘veneer’ of oxidized 
chondritic material which did not segregate into the core 6,9 . In 
each of these hypotheses, metal is added to the surface of the 
growing Earth, and initial interaction between metal and silicate 
will occur at relatively low pressures. Hence, our partition 
coefficients, which are obtained at low pressure, may be appli- 
cable. Later, we shall return to this point and further discuss 
the importance of pressure on partitioning equilibria. 

A fourth hypothesis, that siderophile and chalcophile element 
abundances in the Earth’s upper mantle were established by 
partitioning between metal and silicate at very high pressures, 
perhaps at the core-mantle boundary, will not be explored for 
three reasons. First, metal-silicate partition coefficients have not 
been determined at megabar pressures, precluding quantitative 
evaluation. Second, evidence discussed below suggests that 
siderophile and chalcophile element concentrations in the 
mantle were established no later than 3.5 Gyr ago. Thus, material 
must have been efficiently cycled from the core-mantle boundary 
to the upper mantle and have been well mixed in <1 Gyr. Such 
effective mixing appears problematical. Third, it seems unlikely 
that equilibration of siderophile elements between the core and 
lower mantle would have preserved chondritic ratios of elements 
such as Ni and Co, which differ greatly in their siderophile 
tendencies. 
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Table 1 Experimental conditions and results 


Element 

r(°o 

-•og/o,* 

P(sm/lm> (range) 

£*(ls/lm> (range) 

%s LM t 

N§ 

W 

1,250-1,270 

12.3-12.5 

20(15-30) 

-1(0.12-1.9) 

19-24 

3 

Re 

1,250-1,260 

13.0-13.5 

51 (5 1 - > 60) 

5 x I0~ 4 (5 x t0~ J -5x 10' 5 ) 

23-24 

2 

Ir 

1,270 

12.5-13.2 

46(41-52) 

5 x |0~ 5 (5 x 10~*-5 x 10'*) 

22-23 

2 

Mo 

1,260 



2.45 

— 

25 

0 

Ni 

1,270 

12.5 

1.24 

2x I0“*(4x 10~ 4 -1 x 10 - *) 

23 

II 

Co 

1,260 

12.6-12.7 

2.3 (2.2-2.5) 

0.0075(0.006-0.010) 

24-25 

2 

Au 

1,270 

12.2-13.0 

1.0(0.97-1.1) 

10~ 4 (3x 10~ 4 -5 x 10~ 5 ) 

20-22 

2 

P 

1,250 

12.6 

0.28 

0.24 

18 

l 

Ga 

1,270 

12.5 

3.5 (2.6-43) 

0.8 (0.6-1) 

19-22 

2 

Ag 

1,250 

12.7-13.0 

-0.01 

0.01 (0.01-0.014) 

23-27 

2 

Pb 

1,270 

12.2 

BDLH 

0.15 

20 

1 


* Calculated from the Fe content of the solid metal and the FeO content of the silicate glass using the iron-wustite data of ref. 41. 
t Weight per cent sulphur in liquid metal phase. 

§ Number of individual experiments which contained silicate melt. . .. . . . . ,, 

|| All experiments contain Ni but only one sample (the most oxidized sample to receive a high-neutron-flux irradiation) contained measurable 

Ni in the silicate glass. The range for Z> (LS/LM , is that for two analyses of the same glass, 
f Below detection limit. 


The first three hypotheses will be explored in more detail. 
First, however, we will discuss constraints obtained from experi- 
mental and natural samples which must be satisfied by any 
successful model of core formation in the Earth. 

Constraints 

Partition coefficients. Our experimental partitioning techniques 
have been described elsewhere 10 and will not be repeated in 
detail here. In general, samples containing Fe-Ni metal, natural 
pyrite (FeS 2 ), synthetic basaltic glass and a tracer element are 
placed in an alumina crucible and sealed in an evacuated silica 
tube with a separate Fe-metal-silica assemblage designed to 
maintain the oxygen fugacity of the system (/ 0j , equivalent to 
the partial pressure of oxygen) near the quartz-fayalite-iron 
(QFI) buffer (thus ensuring that iron metal is a stable phase). 
Quenched samples of coexisting solid metal (SM), liquid metal 
(LM) and liquid silicate (LS) are analysed by electron micro- 
probe. If the tracer concentration in the silicate glass is too low 
for microprobe analysis, aliquots of glass are separated from 
metal and analysed by instrumental neutron activation. Table 
1 gives experimental conditions and results for W, Re, Ir, Mo, 
Ni, Co, Au, P, Ga, Ag and Pb— a suite of elements which covers 
a wide range of nebular condensation temperatures (volatilities) 
and geochemical behaviour. In this suite W, Re and Ir are the 
most refractory and Ag and Pb are the most volatile. 
Siderophile abundances in the bulk Earth. The presumed extrac- 
tion of most of the siderophile and chalcophile elements to the 
core makes calculation of the bulk-Earth concentrations of these 
elements difficult or impossible. However, on the basis of 
lithophile elements, which do not readily combine with Fe or 
S, we infer that the bulk Earth is approximately chondritic in 
composition. We do not mean to imply by this that the Earth 
is compositionally identical to the group Cl chondrites, but that 
refractory siderophile elements should be present in chondritic 
relative proportions and that more volatile elements may be 
depleted to varying degrees. Figure l shows the concentrations 
of siderophile and chalcophile elements in several types of 
chondritic meteorites. Refractory siderophile elements are 
usually present in higher concentrations than in group Cl, while 
more volatile elements are typically depleted. Again, although 
the concentrations of siderophile and chalcophile elements in 
the bulk Earth are unknown, it is reasonable to require that 
models of bulk-Earth abundances of these elements be con- 
sistent with siderophile and chalcophile element concentration 
patterns in chondrites. 

Siderophile abundances in the upper mantle. During the past 
decade, a large quantity of high-quality data on mantle samples 
(such as Iherzolites) and basalts has become available. Morgan 




ef a/. 6 and Jagoutz et aL u have presented convincing evidence 
that their Iherzolite samples (or at least their most ‘fertile" 
lherzolite samples) have not undergone large amounts of partial 
melting, and hence are most likely to record primitive mantle 
abundances, at least for the compatible elements. Thus, when 
possible, we have used the mantle abundance estimates of refs 
6 and 11. For the incompatiblle elements P, Mo and W, abun- 
dances are inferred from the trace-element systematics of basalts 
(see below). These estimates are listed in Table 2. 
Mineral-melt partitioning in basalts. The partitioning of Ni and 
Co and between minerals and silicate melts is moderately well 
known from laboratory experiments, but these elements are the 
exception rather than the rule. For most of the elements in our 
suite, solid silicate (SS)/liquid silicate (LS) partition coefficients 
for element i (‘Dss/ls) must be evaluated from element correla- 
tions in natural basalts. For example, a comparison of data for 
basalts and ultramafic xenoliths on Re versus Ir and Au versus 
Ir diagrams from ref. 12 implies that Au ^ss/ls anc * r *£>ss/ls arc 
of order unity, but that U D SS/IJS is —50. (Note in Table 1 that 
during core formation, trivially small fractions of Au, Re or Ir 
enter the silicate phase, making the exact value of these par- 
ticular silicate partition coefficients of little consequence. Also, 
although the host phase of Ir during silicate partial melting is 
unknown, the large differences in partitioning behaviour 
between Au and Ir strongly imply that the host is an oxide or 
silicate.) 

Drake 13 has used element i versus La diagrams to make 
geochemical distinctions between terrestrial, lunar and eucritic 
basalts. Many siderophile and chalcophile elements, when plot- 


Table 2 Siderophile and chalcophile element abundances in the upper 

mantle 


Element 

Concentration/CI* 

Ref. 

W 

0.11 

15 

Re 

0.0075 

6 

Ir 

0.0075 

6 

Mo 

0.064 

15 

Ni 

0.19 

11 

Co 

0.21 

11 

Au 

0.02 

6 

P 

0.058 

14 

Ga 

0.3 

U 

Ag 

0.087 

6 

Pb 

0.09t 

42 


* Concentrations normalized to Cl chondrite ref. 43. 
t Mean of total Pb and unleachable Pb. 
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W Re Ir Mo Ni Co Au P Go Ag $ 204 Pb 


Fig. 1 Sidcrophile and chalcophile element abundances for four 
different types of chondritic meteorites (A, EH; A, EL; •, CV; 
O, H), normalized to the concentrations of these elements in Cl 
chondrites. Elements on the left are more refractory and elements 
on the right are more volatile. Relative to Cl chondrites, refractory 
siderophile and chalcophile elements are typically enriched in 
other chrondrite types and volatile siderophile and chalcophile 
elements are typically depleted. Chondritic materials provide esti- 
mates of the siderophile and chalcophile element abundances that 
may be expected in the bulk Earth (see text). Data are from the 
following sources: E, chondrites, refs 45 (Re, Ir, Ni, Au, Ag), 43, 
46 (Co, Ga), 43, 47 (P), 43, 48 (S), 43, 49 ( 204 Pb); CV, chondrites, 
refs 50 (Ir, Co, Ni, Au, Ga), 43, 51 (Re, Ag), 43, 52 (W), 53 (Mo), 
43, 47 (P), 43, 54 (S), 43, 55 ( 204 Pb); H chondrites, refs 43, 56 (W, 
Re, Ir, Mo, Ni, Co, Au, P, Ga, S) 43, 57 (Ag), 43, 58 ( 204 Pb). When 
only one reference is given, both the chondrite group analysis and 
the Cl normalization analysis are reported in that reference. When 
two references are given, the first is the source of the Cl 
normalization. 


ted on such a diagram, form linear trends which yield informa- 
tion on the depletion of element i in the mantle (relative to 
chondritic abundances) and permit the partition coefficient of 
i between basalt and mantle residuum to be estimated. For 
example, Ga concentrations in basalts increase by a factor of 
two, while La concentrations increase by two orders of magni- 
tude in natural terrestrial samples 13 . If we assume that this basalt 
trend is formed by differing degrees of equilibrium partial melt- 
ing of similar mantle sources and that L "D SS/LS < 10~ 2 , then 
^Dss/ls^OA Correlations of Pb with La show almost exactly 
the same variation as the Ga-La diagram, implying that Pb £>ss/Ls 
is also —0.4. The good (approximately 1:1) correlation of P, W 
and Mo with La or another lithophile incompatible element in 
basalts 14,15 implies that these elements are all very incompatible 
in silicate systems. Diagrams of phosphorous, W and Mo versus 
incompatibile elements are also the best means of ascertaining 
the mantle abundances of these elements 14,15 . Analyses of Ag 
in terrestrial basalts are very scarce; however, based on the 
available data 16,17 , Ag appears to be weakly incompatible like 
Ga and Pb. Consequently, we have adopted a value of A *D SS/LS = 
0.4, although it could be larger. For example, the MORB (mid- 
ocean-ridge basalt) data of ref. 17 are consistent with values of 
A8 £*ss/ls between 0.1 and 0.5, for 5% equilibrium partial melting. 


Summary. Table 3 lists our adopted partition coefficients. As we 
shall see below, the uncertainty associated with these partition 
coefficients, while often large, is of the same order as the uncer- 
tainties regarding the physical conditions under which core 
formation occurred. In this sense at least, we feel that zeroth- 
order modelling, using the data of Tables 2 and 3, is justified. 
We will now examine hypotheses for core formation in the light 
of these data. 

Hypotheses 

Inefficient core formation. The segregation of the core from the 
mantle is envisaged as an equilibration between four phases, 
namely, solid silicate, liquid silicate, solid metal and liquid 
metal, followed by incomplete separation of metal from silicate. 
Equilibrium between these four phases is described by three 
sets of partition coefficients ( D ) for each element at specified 
temperature, pressure and oxygen fugacity: D(solid metal/sul- 
phur-bearing metallic liquid), D(liquid silicate/sulphur-bearing 


Table 3 Adopted partition coefficients 


Element 

®(LS/Lm) 

^(SM/LM) 

Ass/LS) 

W 

1 

36 

0.01 

Re 

5x 10~* 

83 

1 

Ir 

5x 10 -5 

83 

50 

Mo 

8 x 10~ 4 * 

2.45 

0.01 

Ni 

2 x 10 -4 

1.33 

10t 

Co 

7 x 10 -3 

2.3 

3t 

Au 

1 x 10 -4 

1.3 

1 

P 

0.24 

1.7 

0.02 

Ga 

0.8 

6.0 

0.4 

Ag 

0.01 

0.01 

0.4 

Pb 

0.15 

0 

0.4 


* Calculated from our data and that of ref. 25. 

t Ni and Co D iSS/hS) values estimated from ref. 44; all others estimated 
from natural basalt systematics (see text). 

metallic liquid) and D($olid silicate/liquid silicate) where 
is defined to be the weight concentration of element / 
in phase a divided by the weight concentration of element i in 
phase f3. 

The partitioning behaviour of the elements in the experimental 
suite can be generalized in the case of this hypothesis. Mantle 
concentrations of Au, Ir, Mo, Ni and Re are determined almost 
solely by the amount of trapped metal phases. Silver is concen- 
trated nearly exclusively in the trapped metallic liquid; P, W, 
Ga, Co and Pb are intermediate in their siderophilic tendencies. 
Nickel, Co and Ir and compatible (D ss/hs > 1) in solid silicates 
(or oxides); W, P and Mo behave as very incompatible elements 
(D SS /ls« 1); and other elements are intermediate in their sili- 
cate-system compatibility. 

Physical conditions of core formation . Hie fractions of metal and 
silicate in the Earth are taken to be 0.3 and 0.7, respectively, 
the relative masses of the present core and mantle. The propor- 
tions of solid metal and liquid metal are determined by assuming 
that the principal light element in the core is sulphur. Our 
experiments (see below) were conducted so that our metallic 
liquids contained —25 wt% S. Thus, if the present core contains 
8-12 wt % S (ref. 18), use of our experimental data in model 
calculations requires that the metal assemblage which segregated 
to form the core was 30-50% molten at the time of metal-silicate 
equilibration. 

The degree of partial melting of silicate materials during core 
formation is virtually unknown. If the Earth's accretion and 
core formation occurred rapidly, it seems inconceivable that 
high degrees of partial melting could have been avoided 1 ; 
however, the mantle samples themselves give no evidence for 
extensive melting early in the Earth’s history. Spinel lherzolites 
throughout the world are very similar in composition and 
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texture 19 . If Iherzolites arc undepleted in elements that are 
strongly concentrated in basalts, these Iherzolites are described 
as ‘fertile’ (that is, able to produce basalt). Large chemical 
differences arc not observed between the fertile spinel Iherzolites 
and the fertile garnet Iherzolites 19 . Small differences between 
garnet and spinel Iherzolites may exist 20 , but the overall chemical 
similarity between these rock-types implies a mantle which is 
chemically homogeneous down to at least the depth of the 
sources of garnet Iherzolites 150 km; ref. 21). If large degrees 
of partial melting (>20%) occurred in the early Earth, even 
minor density differences between solid and melt would be 
expected to lead to phase separation and the formation of a 
magma ocean 22 . Subsequent crystallization of this hypothetical 
magma ocean would be unlikely to result in a chemically 
homogeneous upper mantle. Rather, it would be expected that 
the formerly molten portion of the mantle would resemble large 
mafic intrusions with pronounced vertical zoning, as observed 
in the Stillwater and Skaergaard intrusions and inferred for the 
Moon 23,24 . Even if such zoning was largely homogenized by 
subsequent mantle convection, it is surprising that apparently 
no remnants remain within the continental lithospheres. The 
uniform chemical composition of mantle Iherzolites (the samples 
from which our mantle siderophile and chalcophile element 
abundances arc mainly derived), coupled with plausible physical 
theories of core formation which predict continuous core forma- 
tion as the Earth grows from the mass of Mars to its present 
mass 3 , lead us to opt for low degrees (<20% ) of silicate partial 
melting. 

The ambient oxygen fugacity during core formation was most 
likely below the iron-wiistite (IW) oxygen buffer (see later 
discussion), but is otherwise imprecisely constrained. Thus we 
allow for isothermal variation in oxygen fugacity in our core- 
formation model. In practice, this means that we allow Dls /L m 
to change with oxygen fugacity in an ideal manner, as defined 
by Rammensee 25 . The *Dls /LM values in Table 3 are appropriate 
for a log fo 2 of -12.75 at 1,250-1,270 °C (approximately one 
log unit below the QFl oxygen buffer). 

Model calculations. The inefficient core formation hypothesis 
may be tested mathematically using the formalism outlined 
below: 

(1) After equilibration among the four relevant phases, most 
of the solid and liquid metal drains away and does not communi- 
cate further with the overlying mantle. This process is subject 
to the constraint of mass balance, so that 

^Earth ** ^SM-^SM + ^LM^LM 

+ f C M X»+ , C u ^L» (!) 

where Xs are mass fractions, l C s are mass concentrations of 
element i and SM, LM, SS and LS refer to solid metal, liquid 
metal, solid silicate and liquid silicate, respectively. Reiterating 


the constraints outlined above, 

Xss + X^ = 0.7 (2) 

+ (3) 

* sm + X lm = 0.3 <«> 

0.3 * X LM /{X LM + X SM ) ** 0-5 (5) 


where p is the fraction of the silicate portion of the system which 
is molten. In terms of partition coefficients (D) and degree of 
partial melting, equation (1) becomes 

= ‘C LM [ i D SM/LM (0.3 - X LM ) + X LM 
+ 0.7 , D ss/lm ( 1 - p) + 0.7 f DLs /LM p] 

As partition coefficients have been measured, *C Earth /*C LM may 
be calculated for assumed values of p and X LM . The sensitivity 
of this calculation to the assumptions of equations (2)-(5) 
should, in general, be acceptably small. For example, the permis- 
sible variation in the total masses of the various reservoirs (solid 
metal, liquid metal, and so forth) is probably about a factor of 


two, and the effects of factor-of-two uncertainties translate into 
two-fold changes in calculated concentrations. 

(2) A small amount of trapped solid metal ( X$ M ) and trapped 
liquid metal (Xl M ) remains behind in the upper mantle. The , 
concentration of element i in the mantle is then given by j 

'Cmtmlc = ^LM^LM + ^SM^SM + ^ss( 1 “ p) + sP (7) 

if X$ M and XJ_ M are small (—0.01). 

Rearranging as in equation (6), we obtain 

‘Cman'to = ^LM (^Lm + ^SM/LM^SM ^ 

+ 1 &ss/lm( 1 — p) + l D LS/LM P) I 

(3) As *C mimU is known from mantle nodules and basalts, 'C LM f 
may be calculated for assumed values of p, X' LM and Xg M . If 
'C LM is known then 'C Ear1h can also bcxalculated from equation \ 
(6). The credibility of the chosen values for the model parameters ’ 
may be evaluated by comparing the calculated to the ; 

concentration of element i in Cl chondrites CCa). Values of | 
f C Eart k for refractory elements which deviate by an order of f 
magnitude from 1-3 x are defined as unacceptable, and new | 
values for the model’s parameters are chosen. > 

There are five important adjustable parameters: X LM , p, / 0j , | 
Xlm and X$ M . However, it must be emphasized that it is 
unknown in detail if the measured partition coefficient values * 
are exactly appropriate. We note that low-pressure ( 1 bar) parti- l 
tion coefficients may be relevant, in that new metal accreting to 
the Earth is deposited at the surface, and metal-silicate equili- , 
brium is likely to be achieved at relatively shallow depths before ■ 
the metal largely segregates into the centre of the growing planet. ~ 
Nevertheless, we should probably be content if any set of model £ 
parameters can be found which, assuming an Earth with chon- j 
dritic relative proportions of refractory elements (2-3 x Cl), will [ 
predict mantle abundances of siderophile and chalcophile ele- j- 
ments to within a factor of two. Better fits, while mathematically : 
pleasing, may have no physical significance. * 

Results. Jones and Drake 7 have evaluated inefficient core forma- £ 
tion as a means of reconciling the high abundances of siderophile l, 
and chalcophile elements in mantle materials with a reducing 
(that is, Iow-/ Oa ) core-forming event. Early results were 
encouraging, but increasing the number of constraints on the ~ 
model (that is, increasing the number of elements to be »J 
modelled) has led to progressively poorer model results. Figure } 
2 shows a ‘best fit’ calculation of bulk-Earth concentrations for p 
our suite of siderophile and chalcophile elements. The refractory ' 
elements W to Mo have concentrations of ~2xCI in the bulk y 
Earth, a reasonable value. The more volatile elements show \ 
various degrees of depletion. Au and P concentrations appear 
low compared with those of more volatile elements. ^ 

Calculated Co and Ni abundances are such that the Co/Ni y 
ratio in the Earth is within 10% of chondritic (a requirement 
of our model), but this agreement is somewhat contrived as Co r 
and Ni are not concentrated in the same phase. Nickel is 
contained mainly in the metallic phases, whereas Co is retained j j 
by the silicates. Maintaining a chondritic Ni/Co ratio in the f- 
bulk Earth by our model requires trapping large amounts of * 
metallic liquid (—2.5 wt % ) in the mantle. Although this amount r. 
of trapped liquid is not inconsistent with theoretical estimates 26 , j 
the retention of this S-bearing liquid leads to unreasonably high y 
S concentrations (~-6,000 p.p.m.) in the mantle, in contrast to i 
the observed concentration of <300 p.p.m. S (ref. 19).'* t 

Furthermore, because we allow for isothermal variation in [ 
oxygen fugacity (and, therefore, in l Dus /LM as well), our best-fit [ 
oxygen fugacity is —1-1.3 log units more oxidizing than the > 
oxygen fugacity necessary for the present mantle FeO concentra- : 
tions to have been estabished by equilibration with metal at J 
I bar. This disagreement is not necessarily a problem, as an j 
oxygen buffer such as QFI will become slightly more oxidizing , 
with pressure (by *-1 log unit per 50 kbar), but more work is , 
needed to define the relative effects of pressure and oxygen r 
fugacity on siderophile and chalcophile element partitioning. ^ 


REVIEW ARTICLE 


NATURE VOL. 322 17 JULY 1986 


W Re I r Mo Ni Co Au P Go Aq S Pb 

Fig. 2 Calculated bulk-Earth abundances of a suite of 12 
sjderophile and/or chalcophilc elements. Elements are listed in 
the same order as in Fig. 1. The model of calculation is that of 
inefficient core formation (sec text). The physical conditions of 
core formation and the amount of trapped metal phases are varied 
until the calculated refractory siderophile element abundances are 
2-3 x Cl for the bulk Earth and the Ni/Co ratio of the bulk Earth 
is within 15% of chondritic. Clearly, geologically reasonable condi- 
tions can be chosen so that these two constraints are satisfied. A 
consequence of the calculation is that volatile elements are found 
to be depleted in the bulk Earth. A drawback of this type of mode! 
is that calculated S abundances for the upper mantle are >20 times 
higher than is actually observed. Best-fit model parameters: trapped 
metallic liquid (25 wt% S), 2.5%; trapped metal, 0.04%; degree 
of silicate partial melting, 10%; log /oj =* - 12.35; S concentration 
in core, 10%. 


Thus, even though it is possible to calculate bulk-Earth con- 
centrations of many siderophile and chalcophile elements which 
are not unreasonable, the consequences of inefficient core forma- 
tion are not totally attractive. If this model is to be viable, there 
must be substantial changes in the values of partition coefficients 
due to an uninvestigated intensive variable (such as pressure). 
Such variations are certainly not out of the question, especially 
as the Fe-FeS eutectic composition is known to be pressure- 
sensitive 27 . 

Equilibrium between an Fe S O metallic liquid and mantle sili- 
cates* This model has been proposed by Brett 8 , and accounts 
successfully for the upper-mantle abundances of a subset of 
siderophile and chalcophile elements. The model postulates 
equilibrium between metal and silicate phases and does not 
resort to trapping small fractions of metal in the mantle. The 
main difference between this model and previous equilibrium 
core-formation models is that no solid metal is present. It recog- 
nizes that liquid metal/silicate partition coefficients may be 
lower than solid metal/silicate partition coefficients (Table 3). 
The redox conditions postulated by Brett* are rather oxidizing, 
in accordance with the absence of solid metal. This model is 
mathematically testable using partition coefficients and mass 
balance constraints in a similar manner to that outlined above. 
Generalizations concerning element partitioning are also as 
above. 

Although at first sight this hypothesis seems successful, there 
are significant difficulties which appear to invalidate it. For 
example, the refractory element Mo was calculated by Brett* to 
be present in the Earth’s mantle at 59 x Cl abundance. Although 
this value is much closer than typical fits obtainable with solid 
metal/silicate partition coefficients, we consider this to be an 
unacceptably high value. In addition, although unknown to 
Brett* at the time, the value of 3.65 x 10” 3 which he used for the 
Ni liquid silicate/ liq uid metal partition coefficient is very large 


(compare with our value of 2x 10“ 4 in Table 3) and, in turn, 
implies an oxygen fugacity too high to be consistent with the 
low FeO content of the Earth’s upper mantle. For example, at 
1 ,250 °C, we predict that the oxygen f ugacity required for Brett’s 
Ni partition coefficient is — 1 log unit more oxidizing than the 
Fe-FeO (IW) oxygen buffer. Under these conditions, the activity 
of Fe in an Fe-S-O metallic liquid will be much less than the 
activity of FeO in any coexisting silicates. The very low iron 
content of the upper mantle (< 10 wt % ) implies that any metallic 
liquid postulated to have been in equilibrium with the upper 
mantle (under oxidizing conditions) had an even lower iron 
activity than the present mantle. 

The possibility that core formation in the upper mantle 
occurred by separation of a metallic liquid with a very low iron 
activity deserves further exploration. In the Fe-Si system at 
1,600 °C, for example, the iron activity at the midpoint of the 
binary system is -0.05 (ref. 28). In the Fe-S system at 1,600 °C 
the iron activity at the midpoint of the system is -0.13 (ref. 29). 
Non-ideal interactions between iron and non-metals can drasti- 
cally reduce the activity of iron, even at fairly high iron con- 
centrations in the liquid. Conversely, decreasing the temperature 
may minimize the non-ideal effects of non-metals. For example, 
at the eutectics of the Fe-S and Fe-Si systems, where the mole 
fraction of iron is definitely less than unity, iron metal is stable 
and the activity of iron in the liquid must also be unity. The 
complex interaction of intensive variables (pressure, tem- 
perature, non-metal concentrations, oxygen fugacity) make 
quantitative prediction of liquid activities difficult. 

Although we cannot rule out highly oxidizing conditions 
during core formation on physical-chemical grounds, we note 
that fractionations of the noble refractory siderophile elements 
are much more probable in this model than in other models of 
core formation. Because noble siderophile elements appear to 
be present in chondritic relative abundances in the Earth’s 
mantle, this is an important consideration. In the inefficient core 
formation model, mantle abundances of the noble refractory 
siderophile elements are ‘buffered’ towards chondritic ratios by 
trapped metal, contrary to the assertion of Morgan 30 . In *late- 
stage veneer’ models , of course, no noble siderophile element 
fractionations are possible because the material is simply added 
to the upper mantle without chemical processing. 

Figure 3 shows the results of recalculating bulk-Earth abun- 
dances for Brett’s* model by extrapolating our internally con- 
sistent set of partition coefficients to higher / 0? . In this model 
AT lm » 0.3, the degree of silicate partial melting is 20%, and f 0j 
is allowed to vary between IW and QFM (quartz-fayalite- 
magnetite, an oxidizing buffer at which iron metal is unstable). 
Our inefficient core formation model results from Fig. 2 are 
shown for comparison. The refractory siderophiles (W, Re, Ir, 
Mo), rather than being enriched relative to Cl, are severely 
depleted. We conclude that Brett’s* equilibrium core formation 
model is no more viable than our own, although this conclusion 
is subject to the same caveat as before — that high-pressure 
partitioning data are virtually non-existent. 

Heterogeneous accretion/‘choDdritic veneer’. The most recent 
advocates of this type of model are Morgan et at 6 and Wanke 9 , 
who have suggested that material accreted to the Earth became 
progressively more oxidized as accretion proceeded. Early 
accretion and core formation was dominated by highly reduced 
materials, and metal effectively segregated siderophile and chal- 
cophile elements to the centre of the growing planet. By the end 
of accretion, an influx of more oxidized materials had caused 
core formation to cease because metal was no longer stable. In 
this model, moderately siderophile elements such as Ni and Co 
stopped segregating to the core when the Earth’s accretion was 
80-90% complete; the noble siderophile elements such as Ir 
and Au ended their segregation into the core when accretion 
was -99% complete. This model predicts that moderately 
siderophile elements should be present in the mantle in chon- 
dritic ratios, at abundances of 0. 1-0.2 x Cl, and that the noble 
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FI*. 3 Calculated bulk-Earth concentrations of siderophile and 
chalcophile elements assuming equilibrium between an Fe-S-0 
liquid and a mantle which is 20% partially molten ( T = 1,250 e C). 
The shaded region shows the calculated concentrations for a range 
of oxygen fugacities, ranging from 1W (upper bound) to QFM 
(lower bound). Bulk-Earth abundaccs calculated by the inefficient 
core formation model are shown for comparison (A). The assump- 
tion that the core formed in relatively oxidizing conditions, without 
any trapping of metals by the mantle, leads to the conclusion that 
the Earth is depleted in refractory siderophile elements, contrary 
to the expectations based on Fig. I. 


siderophile elements should be present in chondritic ratios at 
— 0.01 x Cl abundances. 

This class of hypothesis has gained favour in the past decade, 
in spite of the requirement that extensive physical mixing of 
material into the upper mantle must occur after core formation. 
Impact processes alone seem incapable of such fine-scale mixing 
(witness the heterogeneity of meteorite-derived siderophile ele- 
ments in lunar samples), but sustained sub-solidus mantle con- 
vection may be adequate, particularly if mantle convection is 
layered. The mixing process must be capable of homogenizing 
Fe, Co and Ni at the scale of basalt generation by 3,500 Myr 
bp and be capable of homogenizing Ir at the hand-sample scale 
by today 7,31 . 

In this regard, it is possible that samples of the upper mantle 
which have been stored in continental lithospheres for consider- 
able times (without the benefit of sub-solidus convective mixing 
over geological time) place even tighter constraints on the rate 
of mixing required to homogenize mantle siderophile elements. 
Jagoutz et ai yl have presented evidence that 2,700 Myr-old 
eclogites have been stored in the lower continental lithosphere 
for long periods of time (100-2,000 Myr) without participating 
in mantle mixing processes. More dramatically, Boyd et aL 21 
have argued that eclogitic and garnet peridotitic material, trans- 
ported to the surface by kimberlite eruption, was emplaced in 
the lower continental lithosphere as long ago as 3,000- 
3,500 Myr. Thus, it is quite possible that the isolation of con- 


Table 4 

Test of the multi-component accretion model 


Concentration before 



second episode of 

Concentration after 

Element 

core formation 

removal of metal* 


(C/CI) 

(C/CI) 

Ni 

0.2 

0.2 

Co 

0.2 

0.2 

P 

0.2 

0.2 

Mo 

0.2 

0.064 

Ga 

0.2 

0.2 

W 

0.2 

0.2 


• Just enough metal is removed to lower the Mo concentration from 
0.2 to 0.064 x Cl (see text). 


tinental shield lithospheres from mantle convection and the 
homogenization of Ir that we observe today in continental 
lithospheric materials were both complete by 3,500 Myr ago. 
We question whether mantle convection was capable of such 
rapid, fine-scale mixing. 

Another way of making the same observation is to note that 
Ir contents of Iherzolites from beneath old cratons (such as 
Southern Africa), from beneath younger continental margins 
(such as the southwestern United States) and from mantle over- 
lying oceanic hot-spots (Hawaii) are essentially identical 19 . The 
lithospheres beneath southern Africa and the southwestern 
United States were not obviously isolated from mantle convec- 
tion at the same time, and the lithosphere beneath Hawaii has 
presumably never been isolated to the same degree as continental 
lithospheres. Iridium homogenization was apparently complete 
by the time of the formation of the earliest continental litho- 
spheres. j 

Regardless of absolute rates of mantle mixing, the moderately ! 
siderophile elements are a sensitive test of the heterogeneous j 
accretion hypothesis. Table 4 shows our model results using the * 
partition coefficients of Table 3. For simplicity we have quant- \ 
ized accretion and core formation (before the accretion of the | 
last 1% — the hypothesized 4 Iate-stage veneer’) into two parts: \ 

the first 80% of material accreted is highly reduced, the last f 
20% is more oxidized. This approach permits us to model core [ 
formation as two discrete events: the first, which removed all \ 
siderophile and chalcophile elements from the early-accreted - 
80% of the Earth, and a second event which removed only noble ? 
siderophile elements (such as Ir and Au) from the late-accreted i 
20%. A justification of this simple approach is that, because of 
the large difference in the metal/silicate partition coefficients f 
between Co and Ni, continuous core formation during the last ; 
20% of accretion would be unlikely to preserve the approxi- 
mately chondritic Co/Ni ratio of the upper mantle (Table 2). 
Thus, at the beginning of the second model core-forming event, 
most (if not all) siderophile and chalcophile elements were 
present at ^-0.2 x Cl abundances in the upper mantle, as defined 
by present Ni and Co abundances. At this time we allow just 
enough separation of metal to lower the Mo concentration from 
0.2 x Cl to its present value of 0,064 x Cl (Table 2), and we 
observe the effect of this separation on the other siderophile 
elements. As can be seen from Table 4, no moderately 
siderophile element other than Mo is affected, even though W 
and P are clearly depleted relative to Co and Ni in the present 
mantle (Table 2). 


The heterogeneous accretion hypothesis results in additional 
predictions which are inconsistent with observed elemental 
abundances in the Earth’s upper mantle. First, the amount of 
metal which was extracted from the last 20% of accreting 
material must have been trivial, and the high / Qj requires that 
that metal must have been Ni-rich (because of the decreased 
stability of Fe-rich metal). If we assume that as much as 10% 
of the total Ni was removed during the last core-forming event 
(consistent with the approximately chondritic Co/Ni ratio of 
the upper mantle), then only 0.02% metal was removed. It would 
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volume of (P. T, /<>,)- space relevant to core formation ha* been 
investigated. Two obvious areas for new experimental explora- 
tions are those of high pressure and higher oxygen fugac.ty For 
example from our work on iron meteorites , we have gamed 
some insight into the possible role of pressure. Many s'dcrophile 
trace element-non-metal interactions .n metallic bqu'ds appear 
to be explicable in terms of a simple non-metal-avoidance model. 
Siderophile trace element-sulphur interactions are often 
times stronger than Fe-S interactions, winch are known to be 
pressure-sensitive. (Specifically, the composition of the eutectic 
point in the Fe-FeS binary system becomes more Fe-nch as 


pressure increases 27 .) Extrapolating from the Fe-S system, it is 
possible that siderophile element-non-metal interactions will be 
strong functions of pressure. Quantitatively successful models 
of core formation may require accurate and detailed assessments 
of the (P, T) state of the early Earth. 
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neutron activation analyses; D. Stevenson and F. Richter for 
discussions; and R. Brett, D. J. Malvin, J. W. Morgan and 
H. E. Newsom for helpful reviews. 
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Translocation of proteins across biological membranes is 
a kev step in the intracellular sorting and secretion of proteins. 
Although much is known about protein movement across van0 “* 
membrane systems 1 , it is unclear how hydrophilic, charged 
proteins translocate through the hydrophobic core of a phos- 


Mdhotrexate, a folate antagmist,hU)c]^mp^t into m^o^ondna of mouse d^droft^^e^e^ mo/e , y u does not 

sr that dihydrofo,ate 

must at least partly unfold in order to be transported across mitochondrial membranes. 


pholipid bilayer. In particular, it remains unanswered whether 
a protein must unfold during the transport process or whelher 
it can translocate in a folded state t perhaps through 
non-bilayer domains in the target membrane . 
approached this question by fusing the presequence of an 


